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ABSTRACT 



The dependence of the ratios of the cross sections of deep inelastic lepton-nucleus scat- 
tering is studied in the framework of leading twist, lowest order perturbative QCD. The logQ^ 
slope of the ratio F2'^ / F2 is computed by using the DGLAP evolution equations, and shown 
to be sensitive to the nuclear gluon distribution functions. Four different parametrizations for 
the nuclear effects of parton distributions are studied. We show that the NMC data on the 

dependence of F2^ / F2 rule out the case where nuclear shadowing (suppression) of gluons 
at X ~ 0.01 is much larger than the shadowing observed in the ratio F.^/F^. We also show 
that the possible nonlinear correction terms due to gluon fusion in the evolution equations do 
not change this conclusion. Some consequences for computation of RHIC multiplicities, which 
probe the region x>0.01, are also discussed. 



kari . eskola@phys .jyu.fi 
^heli. honkanen@phys.jyu.fi 
^ vesa.kolhinen@phys.jyu.fi 
^ Carlos . salgado@cern . eh 



1 Introduction 



Nuclear parton distributions (nPDF) are needed in the computation of inclusive cross sections 
of hard, factorizable, processes in high energy nuclear collisions. In the framework of collinear 
factorization and leading twist, it is possible to extract universal nuclear parton distributions 
fi/A{x, Q^) from the measurements of deeply inelastic lepton-nucleus scattering (DIS) and hard 
processes in pA collisions such as the Drell-Yan process. The power corrections in the cross 
sections ^ and in the evolution equations 0, ^ can be neglected if the scales and momentum 
fractions x involved are large enough. The dependence of the nPDF on the scale is then 
given by the DGLAP evolution equations 0] . Analogously to the global analyses of the parton 
distributions of the free proton, the nPDF can be determined based on the DGLAP evolution, 
sum rules and fits to the data. Sets of nPDF like EKS98 ^ (the code can be found in |0, H) 
and HKM § (Hirai, Kumano and Miyama, the code in |jlO|) have become available. 

In the DGLAP analyses of parton distributions, the problem boils down to fixing the initial 
distributions for the DGLAP evolution at an initial scale Q^. For the nPDF in particular, there 
are some uncertainties in the initial distributions due to the lack of experimental data. For 



instance, in the EMC region (a:>0.2), gluons and sea quarks are so far not well constrained [IT 



(see also |T^, Q for more discussion). Also, it seems to be a common belief that the nuclear 
gluon distributions are constrained only very weakly by the DIS data at small values of x. In 
this letter our aim is to emphasize that this is not the case but very valuable constraints for the 
nuclear gluon distributions in the shadowing region [x ~ 0.02) can be obtained from the 
dependence of the structure function ratio -Ff "^Z -^2^ measured by the New Muon Collaboration 
(NMC) ||13[. As first discussed in [|1^, this is a consequence of the fact that in the lowest- 
order DGLAP evolution dF2{x,Q'^) / dlogQ"^ ~ asXg{2x,Q'^) [|l5l at small x. One of the tasks 
in the present paper is to compare the results of the HKM analysis 0, where the constraint 
from the measured dependence of -Ff '^Z has not been applied, with the results obtained 
with EKS98 [^, ^], which makes use of this constraint. 

The Relativistic Heavy Ion Collider (RHIC) initiated the collider era for the search of the 
Quark-Gluon Plasma in ultrarelativistic AA collisions. From the point of view of the nPDF, 
the measurements of charged particle multiplicities dNch/ di] in Au+Au collisions ^7\, [1^, [1^ 
have been truly exciting, since production of semihard gluons at 1...2 GeV scales is expected to 
dominate particle and entropy production at ^/s>200 GeV. The measured multiplicities are thus 
a probe of the nuclear gluon distributions at x ~ 0.01 at these scales: in the models employing 
saturation 0, ^ |21|, ^ |23|, the multiplicity scales as dNch/di] ~ xgA{xsat, QLt) (with possible 



powers of added). In the two-component (hard+soft) models, such as HIJING p4|, the 
nPDF are probed through the perturbative (hard) minijet component. 

Recently, in p5|, it has been suggested based on the HIJING model that the RHIC data 
on multiplicities would indicate that the gluons were more strongly shadowed than the sea 
quarks. To challenge this interesting suggestion, we shall perform a DGLAP analysis of the 
nPDF based on the initial nuclear effects for quark and gluon distributions as given in . As 
the second and main point of this paper, we shall show that in the leading twist, lowest order 
DGLAP framework the NMC data on the dependence of / F2 ||13| rules out very strongly 
shadowed gluons. In relation to HIJING, consequences of this observation for the parameter 
P05 which determines the division into soft and hard components, will be discussed. 
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When the density of gluons in the wave functions of the colhding nuclei (or hadrons) becomes 
large enough, gluon fusion starts to play a role. This induces non-linearities into the QCD 
evolution equations ^, ^ ^ at small values of x and . As the last subject to study 
in this letter, we add the non-linear terms (GLRMQ) into the DGLAP equations for Sn and 
C nuclei. Comparison with the NMC data |]l3| shows, that a very strong gluon shadowing is 
ruled out even more clearly when the non-linearities are included. 



2 DGLAP analysis and different parametrizations 

Following the notation in our previous works |^, , we define the nPDF through the nuclear 
effects, 

Rf{x,Q'^) = ^jj-^^J-, i = 9, u, d, s, M, J, s, (1) 



where ff = f[^^ is the number density distribution of a flavour i in a bound proton of a 
nucleus A, and fi is the corresponding distribution in the free proton. The parton distributions 
of bound neutrons in isoscalar nuclei are obtained through isospin symmetry, f^^^ = f^l^^ and 

fd{i^ = fu[u)^ and we expect this to be a good approximation for non-isoscalar nuclei as well. 
Nuclear effects in the structure function we define through 
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where F2 of the protons (neutrons) of the nucleus A is F2^"'^^^{x, Q"^) = J2q ^l[xfq^"'^^^{x, Q"^) + 
xf^^"'^^^{x, Q^)], and where the small nuclear effects in D have been neglected. 

In order to explicitly show the effect of nuclear gluon shadowing to the scale dependence 
of the ratio Fg^'^/Fg^, we shall study different parametrizations for the nuclear effects. On one 
hand, we will directly make use of the available results from the global DGLAP analyses of 
nPDF, EKS98 [§ |[ and HKM ||]. On the other hand, we perform the DGLAP evolution 
for nPDF by taking the initial modifications Rf{x, Ql) from two different parametrizations in 
which the dependence is assumed to be negligible: 



HPC parametrization |29] which is based on 



/I _ 3; \ Pf 

+ iaerac-bcmcX)Q{x-Xsh)Q{Xf-x)+Rfi-^—-) ^i^-Xf) (3) 

where the different regions are matched together by setting i?sh = Qemc — ^emca^sh, Rf = 

Oemc - bemcXf and a^mc = 1 + &emca;emc- The A dependence of 6emc is &emc = Peine [1 - 

_ 1.145^-2/3 ^ o.93A"i + 0.88^-^/3 _ 0.59^-^/3] ^^^^ ^ fit to DIS data 

results in pA = 0.10011, ca = 0.0127343, cd = 1.05570 x^h = 0.154037, Xeme = 0.275097, 
Peme = 0.525080, Xf = 0.742059, and pi = 0.320992. 
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New HIJING parametrization which replaces the obsolete one in HIJING [E^ (see also 



the discussion in []T2[)- In the new parametrization, the quark sector is fitted to modern 
DIS data, while the gluon sector, especially shadowing and antishadowing are constrained 
by the requirement that (with the updated y/s dependence of the cut-off parameter po for 
the transverse momentum of minijets) HIJING reproduces the measured charged-particle 
multiplicities in Au+Au collisions at RHIC. In this way, gluon shadowing is suggested to 
be much stronger than that of quarks: 

Rf{x)=Rf = Rj^{x) = 1.0 + 1.191og^/^A(x^- 1.2x2 + 0.21x) 

-Sg {A^/'^ - - 3.5v^) exp(-xVO.Ol) (4) 

Rf{x) = 1.0 + 1.19 log^^^A (x^ - l.2x^ + 0.21x) 

-Sg {A^/^ - 1)°-^(1 - 1.5x°-3') exp(-xV0.004), (5) 

with Sq = 0.1 and Sg= 0.24. . . 0.28. Below, we shall use the value 0.24. 

In Fig. ID, in the left panels, we plot the nuclear effects for gluons in tin {A = 117, top 
row) and carbon [A = 12, third row) nuclei at a scale Ql = 2.25 GeV^, as given by EKS98 
and HKM, and by the HIJING and HPC parametrizations above. For further discussion, we 
also show Rp^{x,Ql) from Eq. (|^) (2nd and 4th rows, correspondingly). With the EKS98, 
HIJING and HPC parametrizations, we have used the MRST (central gluon) LO PDFs of the 



free proton p2| , ^. The HKM results are obtained directly from the HKM code |T0[, based on 
and where the LO MRST (c-g) PDFs are also used. 

As seen in Fig. |1], for the gluons there are quite distinctive differences between the sets 
used: in EKS98 there is strong gluon antishadowing, which originates from the requirement 
of conservation of momentum combined with the constraint obtained from the measured 
dependence of F2^/F2- In HKM in turn, the gluons are less shadowed at small values of x 
but no EMC-effect appears at a; ~ 0.3 . . . 0.7. Momentum is conserved also in HKM, so the 
deficit of momentum at x<0.2 is compensated by a rapid increase of R^ at x>0.2. Of the four 
cases studied here, the HIJING parametrization has the strongest gluon shadowing. Since no 
antishadowing appears, the HIJING parametrization underestimates the momentum sum at 
by about 10 %. As discussed in |l^, the HPC parametrization underestimates the momentum 



sum at Qq by 5 %. 

Regarding the ratio Rp^-, all parametrizations are quite similar at x>0.02, as they are based 
on the fits to the DIS data. At x<0.005, however, due to the lack of DIS data in the perturbative 
region, some differences arise. Note also that in Fig. |l| the ratio Rp^ from EKS98 is plotted 
for isospin symmetrized A = 117 (in order to compare with the NMC data [|l3l) but the HKM- 
results are for non-isospin symmetric tin, and with the small effects for D included, as obtained 
directly from the numerical code |10| by HKM. At the region of small values of x, where the 



focus of the present paper is, the isospin effects can in any case be safely neglected. 

Next, we consider the (lowest order) DGLAP evolution of the nPDF and explicitly show 
the consequences of different assumptions of gluon shadowing. In EKS98 and HKM the scale 
evolution has already been done but for the other two cases it needs to be performed separately. 
For the HPC and HIJING nuclear effects, we do this by choosing Qq as the initial scale, and 
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Q- = 2.25 GeV 



= 100 GeV^ 




Figure 1: Rg{x,Q'^) and Rf^{x^Q'^) for Sn {A = 117) and C = 12) as functfon of x for scales 
Q'^ = Ql = 2.25 GeV2 and = iqo GeV^. The difference of the ratios Rf{2x,Q^)/R^^{x,Q^) for 
Sn and C is proportional to the slope d{F2^ / F2) / dlogQ"^ ^ as shown in Eq. @. 



computing the intial distributions with the MRST distributions, taking the initial nuclear effects 
from Eqs. (0) and (^-(§), correspondingly. The absolute nPDF are then evolved from Ql to 
higher scales with the DGLAP equations. The results for Rq{x,Q'^) and Rp^{x,Q^) at a scale 

= 100 GeV^ from all cases studied, are shown in the right panels of Fig. |l|. Again, with 
the EKS98 we use the MRST distributions in plotting the ratio Rp^, and the HKM results are 
directly from the HKM code . 

Two observations can be immediately made: first, the pQCD scale evolution is a sizable 
effect for the gluon ratios at a;<0.3, especially in the small-x region. Second, at small values 
of X, due to the very strong gluon shadowing, the HIJING parametrization predicts the slope 
dRp^{x, Q^)/d log Q"^ to be negative, contrary to the other three cases studied. 

The comparison between the calculated results and the NMC data for the ratio -Ff F2 
is shown in Fig. 0. The EKS98 parametrization reproduces well the experimental data, the 
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Figure 2: Comparison of the calculated and measured dependence of the ratio ■ The 

NMC data |13| are shown with statistical errors only. The results for EKS98 j^, ^ (solid lines) and 
HKM Q (dotted-dashed) are from the corresponding global DGLAP analyses. The dependence of 
the HPC (dashed) and HIJING (dotted) cases is obtained from the DGLAP equations by taking the 
initial conditions for the nuclear effects from Eqs. (^) and (0)-(|5|). 



log Q"^ slopes in particular, which is not surprising as these data have been taken into account 
in the analysis ^. The DGLAP-evolved HPC also reproduces the data reasonably well. The 
HKM results miss the absolute normalization of the data at the smallest values of x but have 
the right kind of curvature. The data would also suggest faster evolution at small scales than 
that in HKM. The HIJING parametrization results in a negative Q^-slope for the ratio F^^ jF^ 
at small values of x and . This clearly is in contradiction with the data. This behaviour is 
caused by the strong gluon shadowing in the HIJING parametrization at small x, as will be 
discussed next. 

At small values of x, the structure function F^ (of both p and n) is dominated by sea 
quarks, whose DGLAP evolution in turn is dominated by gluons. The logQ^ slope of F2 can 
be approximated at lowest order as |TB 



dFt-\x,Q^) 

aiogg2 



xg{2x, 



This leads to 



ai?4(x,Q2) l0asXg{2x,Q' 



aiogQ2 27n lFP{x,Q^) 



Rti2x,Q^ 



x,Q^ 



(6) 



(7) 
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and 



aiogg2 2777 lF^ix,Q^)R'^^{x,Q^) \ R%{x,Q^) R%{x,Q^ 

where xg is the gluon distribution in the free proton, and = + FJ'. 

In the HIJING parametrization gluons are much more strongly shadowed than quarks, so 
R^{2x,Ql) < Rg{x,Ql). Eq. (|^) thus directly shows why a negative logQ^ slope for the ratio 
Rp^ seen in Fig. || is bound to follow. Based on Eq. (H) we can also understand the origin 
for the differences between the computed logQ^ slopes of the ratio F2^/F2- The second term 
in the curly brackets in Eq. is obviously always closer to unity than the first one, since 
the nuclear effects in smaller nuclei are smaller. On the other hand, the NMC data in Fig. 
indicates a clearly positive logQ^-slope. These facts imply that R^ is bounded from below as 
R^{2x,Q^) > Rp^i^XjQ"^). Consequently, R^{x,Ql) cannot be much smaller than Rp^{x,Ql). 
The values of R^{2x,Ql) and Rp^{x,Ql) at x = 0.0125 for Sn and C are directly readable off 
from Fig. ^ Both the sign and the relative order in the magnitude of the computed slopes 
in the first panel of Fig. H can be understood by substituting these values in Eq. (|]). The 
fact that the log slopes from the HIJING parametrization of nuclear effects are opposite to 
the measured ones leads us to the conclusion that very strongly shadowed gluons are ruled out 
within the leading twist DGLAP framework. 

The observation above has an interesting consequence for the HIJING model. As discussed 
in I^Sf, particle production in AA collisions in HIJING is due to contributions from soft and 



hard components, 

— 1 - oV^^partsA^soft; + Unary) \^hs.rd) jjj^^-j=— , [yj 

ar] A (Jin [y's) 

where (Ai'parts) ~ 2A is the average number of participants in AA for the centrality selection 
considered, (risoft) = 1-6 is the average multiplicity from soft processes to the rapidity inter- 
val considered, (A^binary) ~ ^^^^ is the average number of binary nucleon-nucleon collisions, 
('^hard) = 2.2 represents particle production from (mini)jet hadronization, and o^^^{^/s) is the 
inelastic nucleon-nucleon cross section. The integrated minijet cross section cr^t'^y/sjPo) is 
computed perturbatively by using the nPDF. The parameter po for the minimum transverse 
momentum of minijets defines the division into soft and hard components. The values of Pq are 
determined based on fits to the inelastic cross sections measured in pp and pp collisions. The 
^/s dependence of po in HIJING has now been updated , and po is an increasing function of 



i/i. Unlike the saturation scale in models with parton saturation p|, |20|, |2T|, |2^, |2^ , the scale 
Po does not depend on A. 

We write afJ{y/s,po) = [Rf{{x), (Q^))]P(Xjct(yi,Po), where (x) ^ 2po/v^ and {Q^) = apl 
with a ~ 2 for the minijet production at central rapidities. The effective power p = 1 ... 2 
describes the net effect of gluon shadowing. From Eq. (|^) we then obtain 



2' parts) ("-soft) 



(10) 



binary / 

Requiring that the model reproduces the measured multiplicity in Au+Au collisions at RHIC, 
the r.h.s. of Eq. (|TUp is a fixed number for each -^i and A, and thus independent of po- 
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As discussed above, within the DGLAP framework the NMC data imphes that should 
be larger than that of the HIJING parametrization in Eq. (j^) (this is due to both weaker 
gluon shadowing and due to scale evolution). Correspondingly, crjet{y/s,Po) should be smaller. 
This in turn implies that po should be larger than in the pp case. We are thus lead to the 
conclusion that the same multiplicities as are currently obtained from the HIJING model with 
an ^-independent po{^/s) and very strong gluon shadowing, can be obtained with weaker gluon 
shadowing by introducing a scale po which is an increasing function of both ^/s and A. 



3 Non-linear effects in evolution 

The DGLAP scale evolution discussed above is linear in the parton densities. At very small 
values of x the density of gluons increases to the extent that contributions from non-linear 
phenomena due to gluon fusion may start to play a role. The first non-linear correction terms 
to pQCD evolution equations have been computed in |^, let us call them "GLRMQ terms". 
One could argue that perhaps these corrections would change the scale evolution in such a way 
that stronger gluon shadowing could be allowed. In order to study this possibility, we next add 
the GLRMQ terms into the DGLAP equations. In connection with the DGLAP analysis of 
nPDF, the non-linear effects have been numerically studied e.g. in |Q (see also p^). 

The valence quark evolution remains unmodified, but for the gluons and sea quarks the 
generic form of the equations reads 

dxgA ( dxQA \ 1 ^ ^2^ 



dlogQ^ VSloggVDGLAP Q2-^^^^(^'^) ^^^^ 



dxqA ( dxqA \ 1 



(12) 



d log 



-TlgggiX^Q') (l3) 



where GLRMQ terms TZggg ~ a'^s Ixidy/y)y^g^A\y^Ql)^ T^ggg ~ "^a^^^A^a;, Qg) and 7^?^ 
as Jx{dy/y){x/y)'yFG{x/y)yguT{y,Q'^)- The detailed form of these terms can be found in | 



Following |H] and |33], for the 2-gluon density we take x'^g^2\x,Q'^) = -^[xgA{x,Q^)]'^, where 
^j:^ and the nuclear radius Ra = 1.12 A^^^— 0.86 A~^^^ are based on the Woods-Saxon parametriza- 
tion of nuclear densities. For the gluon higher-twist term, we assume xgnrix, Ql) = x^g^A (x, Qo) 
and that gnxix, Q"^) > 0. 

We obtain the initial conditions for the actual nPDF as before, from the EKS98, HIJING 
and HPC parametrizations and with the MRST distributions for the free proton. We do not 
make an attempt to include the GLRMQ terms to the HKM analysis. The results for the scale 
evolution of the ratio yy7-^2^°/ j2^2 shown in Fig. || against the NMC data. We observe that 
the effects of the GLRMQ corrections remain fairly modest (as they should, in order to stay 
as corrections) but that they make the positive logQ^-slopes (EKS98, HPC) flatter, and the 
negative slopes (HIJING) even more negative than in the case of pure DGLAP evolution. Our 
conclusion therefore is that the GLRMQ terms do not give support to the strong shadowing of 
gluons, either. 
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Figure 3: The scale dependence of the ratio /F^ calculated using DGLAP evolution with MQ 



corrections, and compared to NMC data |13]. Initial conditions for nuclear effects are taken from 
EKS98 (solid lines), HPC (dashed) and HIJING (dotted) parametrizations. 



The systematics of the change in the logQ^ slopes at the initial scale Qq can also be easily 



Ai 



TP A' _L 
-^2 DGLAP "T 



understood from Eqs. (|T^. The evolution equation for F.f can be written as F. 
-^2GLRMQ5 where the prime stands for (9/9 log Q^, and where the second term contains only 
the GLRMQ corrections for the sea quark evolution from Eq. (0). The net effect of the 
GLRMQ corrections to the log slope of sea quarks is negative, and dominated by the term 



d{Ft/F^) 
d log 



The log slope of the the ratio -Ff '^Z F2 can then be expressed as 



FT 
F^ 



-^2 DGLAP 



-^2 DGLAP 



Ft 



F^ 



+ 



pSn/ 

-^2 GLRMQ 

Ft 



-^2 GLRMQ 



F^ 



(14) 



where the latter term in brackets again contains only the GLRMQ contributions to the log 



slope from Eq. (|T2|). Using F2 



{F^ + F2"), the GLRMQ part becomes 



TpSnl 

-^2 GLRMQ 

Ft 



TPCi 

-^2 GLRMQ 



12V3(i?C) 



117 
"12" 




pSn 



(15) 



Reading the values for Rp^{x, Ql) and Rf{x, Ql) at x ~ 0.01 off from Fig. 1, it is easy to see 
that the contribution from the GLRMQ terms to the logQ^ derivative of the ratio Ft/F2 is 
indeed negative in all the three cases studied. 

We emphasize that for the DGLAP+GLRMQ case we have not attempted to make a global 
analysis of the nPDF, where the initial conditions would be based on fits to various set of data. 
Such an analysis will require also detailed studies of the constraints for the magnitude of the 
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nonlinearities in the case the free proton. We note, however, that according to Fig. ^, in the 
framework of DGLAP+GLRMQ we could expect somewhat less gluon shadowing for QaIx, Ql) 
than in the framework of pure DGLAP. This will also lead to slightly smaller excess of gluons 
at larger values of x (see EKS98 and HKM in Fig. 1). This analysis is left as a future task. 



4 Conclusions 

We have studied the dependence of the ratio / F2 in the pQCD framework of lowest- 
order DGLAP evolution and leading twist. We emphasize that the NMC data on the 
dependence of F2'^/F2 |TB| is the most direct (although indirect) measurement of the nuclear 
gluon distributions that is currently available. These data provide a valuable constraint for 
pinning down the nuclear gluon shadowing in the global DGLAP analyses of nPDFs. 

We have demonstrated the sensitivity of the log slopes of -Ff °/ F2 to the gluon shadowing 
by comparing four different approaches. The NMC data [0 implies that the nuclear effects in 
the gluon distributions (-R^ ), should lie close to those in F2 / F2 at a; ~ 0.01 and ~ 2 GeV^. 
In particular, a very strong gluon shadowing, as suggested e.g. in [^, is ruled out by the NMC 



data, since it leads to a log slope whose sign is opposite to what is measured. Consequences 
for the HIJING model have been discussed in Sec. 2. We also suggest to use the NMC data 
[^] and as further constraints in the DGLAP analysis of HKM . 

The DGLAP approaches 0, ^ and conserve momentum explicitly. The scale indepen- 
dent parametrizations from HIJING |^4|, ^ and HPC do not do that, especially not if 



gluons are very strongly shadowed and no antishadowing appears. We emphasize that the 
strong antishadowing in the gluons of EKS98 results from the fact that the gluon shadowing 
is constrained by the NMC data at x ~ 0.01, and to compensate the loss of momentum there, 
antishadowing is needed. 

We have also studied the effects of the nonlinear GLRMQ terms [0, |^ in the DGLAP 
equations. These terms decrease the slope d{F2'^/ F2) / d log Q^. For the very strongly shadowed 
gluons this results in even more negative slopes than without the GLRMQ corrections. Our 
conclusion therefore is that within the DGLAP and DGLAP-I-GLRMQ frameworks studied we 
cannot find any support from the DIS data for a much stronger gluon shadowing at x ~ 0.01 
and ~ 2 GeV^ than what is observed in the ratio F2 I F2 ■ To resolve the situation at smaller 
values of x, especially in the region relevant for the few-GeV scales at the LHC, more DIS data 
in the perturbative region would be needed. 
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